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THE  ANALYSIS  OF  METAL-FORMING  PROCESSES 
E.  H.  Lee  and  R.  L.  Mallet t 

1 .  Evaluation  of  Stress  and  Deformation  in  Metal-Forming  Processes 

A  complete  stress  and  deformation  analysis  of  a  metal-f orraing  process  is 
necessary  in  order  to  assess  the  onset  of  metal-forming  defects  such  as  the 
initiation  of  internal  or  surface  cracks,  the  generation  of  high  residual 
stresses  or  the  occurrence  of  local  flow  abnormalities.  In  order  to  evalu¬ 
ate  stresses  throughout  the  material  it  is  necessary  to  carry  out  an  elastic- 
plastic  analysis.  Although  in  metal  forming  processes  elastic  strains  are 
usually  much  smaller  than  plastic  strains,  rigid-plastic  theory  is  not  appro¬ 
priate  since  it  cannot  predict  stresses  in  the  rigid  regions.  In  the  actual 
material  these  correspond  to  regions  of  purely  elastic  deformation  and  to 
regions  of  contained  plastic  flow  where  the  total  strains  are  constrained  to 
be  of  elastic  order. 

The  required  finite-deformation  elastic-plastic  analysis  has  only  recently 
become  feasible  due  to  developments  in  continuum  mechanics  for  large  strain, 
related  computational  techniques  and  computer  capacity.  A  finite-element 
computer  program  valid  for  elastic-plastic  analysis  at  finite  strain  was 
developed  ancl  was  first  applied  to  a  problem  involving  plane  strain  extrusion 
through  frictionless  curved  dies.  In  order  to  extend  this  new  capability  to 
encompass  additional  features  encountered  in  practical  metal-forming  processes, 
a  method  of  handling  friction  at  metal-tool  interfaces  was  developed.  The 
required  non-self-adjoint  mathematical  description  led  to  a  finite-element 
formulation  involving  non-symmetric  incremental  stiffness  equations  and  thus 
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to  substantial  modification  to  the  computer  program.  The  results  obtained  for 
plane  strain  and  axisymmetric  chamber  extrusion  through  curved  dies  conclusively 
demonstrated,  as  argued  in  the  project  report  [1]*,  the  feasibility  of  the 
required  large-strain  elastic-plastic  analysis  and  constituted  a  major  step 
toward  the  goal  of  rational  assessment  of  practical  metal  forming  processes. 

Numerical  accuracy  and  computational  efficiency  were  emphasized  throughout 
the  development  of  the  computer  program.  In  order  to  efficiently  handle  materials 
exhibiting  the  common  saturation  of  hardening  at  large  strains,  considerable 
revision  of  the  computational  procedures  was  necessary  because  the  flattening 
of  the  stress-strain  curve  generates  numerical  instability  associated  with  the 
incompressible  nature  of  plastic  flow.  Computer  program  modifications  include 
the  introduction  of  a  partial  stiffness  technique  to  permit  raidstep  yielding, 
of  a  predictor-corrector  procedure  to  establish  midstep  stresses  and  yield  sur¬ 
face  normals,  of  an  exact  procedure  for  tracking  the  strain  hardening  curve, 
and  of  a  constant  dilatation  isoparametric  quadrilateral  finite-element,  as 
well  as  several  higher  order  elements,  to  improve  accuracy  in  the  computed 
stress  distributions.  The  net  result  was  a  marked  improvement  in  accuracy  for 
a  given  spacial  mesh  with  no  appreciable  increase  in  computation  time. 

As  part  of  an  international  effort  to  compare  and  verify  existing  elastic  - 

plastic  computer  codes,  the  upsetting  of  a  cylindrical  billet  between  rigid 

parallel  platens  was  analyzed  and  presented  in  project  report  [2].  We  have 

assumed  responsbility  for  planning  and  coordinating  the  next  phase  of  the 

effort.  We  have  compared  finite-element  and  slip-line  field  theory  solutions 

for  a  related  problem,  plane-strain  compression  of  a  rectangular  block,  in 
_ 

Numbers  in  square  brackets  refer  to  the  bibliography  at  the  end  of  the  report. 


which  large  strains  at  the  block's  edges  produce  folding  of  the  free  surfaces 
against  the  advancing  platens.  This  phenomenon  generates  additional  compu¬ 
tational  difficulties.  In  order  to  analyze  such  problems  satisfactorily 
using  finite  elements,  the  boundary  conditions  had  to  be  modified  to  permit 
boundary  shear  slip  in  the  case  of  perfectly-plastic  material. 

2 .  Incremental  Elastic-Plastic  Theory  at  Finite  Deformat ion 

The  plastic  part  of  an  elastic-plastic  deformation  is  that  remaining  when 
the  stress,  and  hence  the  elastic  strain,  is  reduced  to  zero.  Elastic  defor¬ 
mation  is  that  produced  in  this  un-stressed,  purely-plastically-deformed  ma¬ 
terial  by  the  action  of  stresses  up  to  yield.  The  exact  mathematical  repre¬ 
sentation  of  the  kinematics  of  such  elastic-plastic  deformation  at  finite 
strain  and  rotation  expresses  the  total  deformation-gradient  matrix  as  the 
plastic  deformation-gradient  matrix  pre-multi.plied  by  the  elastic  deformation- 
gradient  matrix  [3].  The  fact  that  this  operation  is  not  commutative  immedi¬ 
ately  throws  doubt  on  the  almost  universal  assumption  in  elastic-plastic  theory 
that  the  total  strain-rate  is  equal  to  the  sum  of  elastic  and  plastic  strain- 
rates,  for  this  relation  is  commutative.  An  exact  analysis  based  on  the  non¬ 
linear  kinematics,  [4],  was  developed  which  shows  that  the  theory  in  current 
use  is  in  error  to  the  extent  of  not  being  objective,  i.e.  deductions  from  it 
may  depend  on  the  particular  axes  chosen  to  express  the  deformation. 

In  order  to  generate  an  analytically  rigorous  theory,  the  nonlinear, 
finite-deformation,  isotropic ,  objective,  elastic  law  is  used  relating  the 
stress  to  the  elastic  deformation  gradient.  In  the  development  of  the  incre¬ 
mental,  or  rate  type,  elastic-plastic  relation ,  the  elastic  law  for  stress  is 


differentiated,  yielding  rate  equations  which  are  therefore  also  objective. 

This  follows  directly  from  the  mathematical  manipulation  thus  avoiding  the 
need  to  select  a  particular  stress-rate  definition  [4].  The  elastic  rate 
equation  and  the  classical  plastic  rate  relation  are  substituted  into  the 
nonlinear,  finite-deformation  kineraatical  theory  to  generate  an  elastic-plas¬ 
tic  constitutive  relation  involving  rate  of  stress,  rate  of  total  strain,  the 
stress  and  the  work-hardening  based  on  deformation  history.  This  constitutive 
relation  exhibits  the  symmetry  properties  to  give  a  rate-potential  function  [5] 
/in'!  hence  a  convenient  variational  principle  for  finite  deformation  analysis 
using  the  finj  te-eierrent  method  [4].  The  theory  is  contrasted  with  that 
in  common  use  and  anomalies  in  the  latter  are  discussed.  The  finite-deformation 
elastic-plastic  theory  in  current  use  is  generated  as  an  approximation  to 
the  rigorous  theory,  and  it  is  anticipated  that  this  approximation  will  be 
adequate  for  most  metal-forming  problems.  However,  it  is  important  from  a 
basic  standpoint  to  know  which  theory  is  rigorous  and  which  is  approximate,  and 
there  may  be  some  problems  of  plastic  instability  for  which  the  difference  is 
significant.  For  situations  involving  high  hydrostatic  pressure,  such  as  occur 
in  shock-waves  and  high-velocity  impact;  where  finite  dilatational  elastic  strains 
arise,  the  new  theory  is  needed  to  provide  an  adequate  analysis  [6,7,8]. 
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A  complete  stress  and  deformation  analysis  of  a  metal-forming  process  is 
necessary  in  order  to  assess  the  onset  of  metal-forming  defects  such  a.  the 
initiation  of  cracks,  the  generation  of  high  residual  stresses  or  the  occur¬ 
rence  of  local  flow  abnormalities.  Rigid-plastic  theory  can  only  predict 
stresses  in  the  regions  exhibiting  significant  current  plastic  flow  so  that 
to  evaluate  stresses  throughout  the  material  it  is  necessary  to  carry  out  an 
elastic-plastic  analysis.  A  finite-element  computer  program  to  evaluate  com¬ 
plete  stress  and  deformation  distributions  has  been  developed  and  applied, 
bringing  new  insights  to  the  assessment  of  metal-forming  processes. 

The  plastic  part  of  an  elastic-plastic  deformation  is  that  remaining  when 
the  stress,  and  hence  the  elastic  strain,  is  reduced  to  zero.  Elastic  defor¬ 
mation  is  that  produced  in  this  purely  plastically  deformed  material  by  the 
action  of  stresses  up  to  yield.  The  associated  exact  finite-deformation  kine¬ 
matics  shows  the  almost  universal  assumption  that  the  total  rate  of  deformation 
is  the  sum  of  elastic  and  plastic  rates  to  be  in  error.  An  incremental  elastic- 
plastic  theory  is  developed  using  the  nonlinear  kinematics.  The  theory  is  con¬ 


trasted  with  that  in  common  use  and  anomalies  in  the  latter  are  discussed. 


